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ABSTRACT 

Observations and theory suggest that star clusters can form in a subvirial (cool) state 
and are highly substructured. Such initial conditions have been proposed to explain 
the level of mass segregation in clusters through dynamics, and have also been suc- 
cessful in explaining the origin of trapezium-like systems. In this paper we investigate, 
using iV-body simulations, whether such a dynamical scenario is consistent with the 
observed binary properties in the Orion Nebula Cluster (ONC). We find that several 
different primordial binary populations are consistent with the overall fraction and 
separation distribution of visual binaries in the ONC (in the range 67 - 670 au), and 
that these binary systems are heavily processed. The substructured, cool-collapse sce- 
nario requires a primordial binary fraction approaching 100 per cent. We find that 
the most important factor in processing the primordial binaries is the initial level of 
substructure; a highly substructured cluster processes up to 20 per cent more systems 
than a less substructured cluster because of localised pockets of high stellar density in 
the substructure. Binaries are processed in the substructure before the cluster reaches 
its densest phase, suggesting that even clusters remaining in virial equilibrium or un- 
dergoing supervirial expansion would dynamically alter their primordial binary popu- 
lation. Therefore even some expanding associations may not preserve their primordial 
binary population. 
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1 INTRODUCTION 

It is thought that the vast majority of stars form in clustered 
environments (with surface densities of several stars, to sev- 
eral h u ndred s tars per square parsec, e.g.lLada fc Lada 2003; 
lLadal l201fj iPortegies Zwart. McMillan fc Gielesl |2010| ). 
Whether all such clusters are dense enough to dynamically 
process the primordia l stellar population is currently the 
subject of debate (e.g. lBressert et al]|2010l ). However, there 
is observational and theoretical evidence that some clusters 
do at least undergo a dense phase in their evolution, a no- 
table example bein g the Orion Nebula C luster (ONC). 

Recent work bv lAllison et al. I (|2009l ) has shown that the 
observed mass segregation in the ONC can be o f a dynamical 
origin . If a cluster is initially substructured (|Allison et al.l 
(2009) used fractals to create substructure) and subvirial, 
then the cluster undergoes cool-collapse and the most mas- 
sive stars mass segregate, in some cases forming trapezium- 
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like systems (| Allison fc Goodwin! liuTlh . Previously, it had 
been thought tha t the mass segregation in the ONC had 
to be primordial (|Bonnell fc Daviea Il998l ). as the level of 
dynamical mass segregation required cannot occur within 
~ 1 Myr in clusters with smooth radial profiles. 

Given the success of the cool-collapse model in pro- 
ducing the observed levels of mass segregation and trapez- 
ium systems, an investigation into the effects of this 
dynamical scenario on clusters containing p rimordial bi- 
nary populations is timely. For simplicity, lAllison et al.l 
(2009) did not include primordial binaries in their simu- 
lations. However, the binary fraction in the O NC is con- 
sistent with that in the field (~ 45 per cent, iPetr et al 
1 199: " 



in the field 

2007 ). A^-body simu l ations bv lKroupa 



Reipurth el al 
( 1995al lb h: iKroupa. Petr fc McCaughreanl (|l999l ). and 



recently, iParker et al.l ((2009), have shown that in a dense 
cluster in virial equilibrium, a binary population with a high 
primordial binary fraction (~ 100 per cent) will be processed 
to a much lower bin ary fraction, c onsist ent with the obser- 
vations in the ONC. IParker et all (2009) argued that it was 
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unlikely that the primordial binary population in the ONC 
was field-like, as the ONC is expanding (indicating that it 
was much denser in the past and therefore had a higher pri- 
mordial binary fraction) and there are no wide (>10 00 au) 
binary systems l|Scalfv. Clarke fc McCaughreanll 19991 ) (indi- 
cating that the binary population has been well processed). 

In this paper, we investigate the effect of dynamical 
evolution in substructured, subvirial clusters on various pri- 
mordial binary populations. We run suites of Af-body simu- 
lations in which we vary the initial amount of substructure, 
and the proportion of stars in binary systems, and compare 
the results to the most recent observations of binaries in the 
ONC. In Section [2] we describe the set-up of the clusters, 
and the initial binary populations; we present our results 
in Section [3] we provide a discussion in Section [4] and we 
conclude in Section [5] 



2 METHOD 

2.1 Initial conditions 

The clusters we simulate have 1500 members, which corre- 
sponds to a cluster mass of ~ 10 3 Mq . For each set of initial 
conditions, we run an ensemble of 10 simulations, identical 
apart from the random number seed used to initialise the 
positions, masses and binary properties. 

Our clusters are set up as fractals; observations of 
young unevolved star forming regions indicate a high 
level of substructure is present (i.e. they do not have 



a radially smooth profile, e.g. [ Cartwrig ht fc Whitworthl 
12004 ISanchez fc Alfaroll201o1 : ISchmeiall201ll . and references 
therein). The fractal distribution provides a way of creating 
substructure on all scales. Note that we are not claiming 
tha t young star clusters are fracta l (although they may be, 
e.g. lElmegreen fc Elmegreenll200ll ). but the fractal distribu- 
tion is a relatively simple method of setting up substruc- 
tured clusters, as the level of substructure is described by 
just one parameter, the fractal dimension, D. In three di- 
mensions, D = 1.6 indicates a highly substructured cluster, 
and D = 3.0 is a roughly uniform sphere. 

We set up the fractals according to the method in 
iGoodwin fc Whitworthl l|2004l ). This begins by defining a 
cube of side N div (we adopt N div = 2.0 throughout), in- 
side of which the fractal is built. A first-generation parent 
is placed at the centre of the cube, which then spawns A"| iv 
subcubes, each containing a first generation child at its cen- 
tre. The fractal is then built by determining which of the 
children themselves become parents, and spawn their own 
offspring. This is determined by the fractal dimension, D, 
where the probability that the child becomes a parent is 
given by A^ -3 '. For a lower fractal dimension fewer chil- 
dren mature and the final distribution contains more sub- 
structure. Any children that do not become parents in a 
given step are removed, along with all of their parents. A 
small amount of noise is then added to the positions of 
the remaining children, preventing the cluster from having a 
gridded appearance and the children become parents of the 
next generation. Each new parent then spawns N^ iv second- 
generation children in AT|; V sub-subcubes, with each second- 
generation child having a probability of becoming 
a second generation parent. This process is repeated un- 
til there are substantially more children than required. The 



children are pruned to produce a sphere from the cube and 
are then randomly removed (so maintaining the fractal di- 
mension) until the required number of children is left. These 
children then become stars in the cluster. 

To determine the velocity structure of the cloud, chil- 
dren inherit their parent's velocity plus a random compo- 
nent that decreases with each generation of the fractal. The 
children of the first generation are given random velocities 
from a Gaussian of mean zero. Each new generation inher- 
its their parent's velocity plus an extra random component 
that becomes smaller with each generation. This results in a 
velocity structure in which nearby stars have similar veloci- 
ties, but distant stars can have very different velocities. The 
velocity of every star is scaled to obtain the desired virial 
ratio of the cluster. 

We set up clusters with fractal dimensions of D = 1.6 
(very clumpy), D = 2.0 and D = 3.0 (a roughly uniform 
sphere), in order to investigate the full parameter space. 
The clusters are out of virial equilibrium at the start of 
the simulations and have a virial ratio of Q — 0.3, where 
we define the virial ratio as Q = T/|fl| (T and |fi| are the 
total kinetic energy and total potential energy of the stars, 
respectively). Therefore, a cluster with Q = 0.5 is in virial 
equilibrium and a cluster with Q = 0.3 is 'subvirial', or 
'cool'. 

To create a stellar syst em, the mass o f the primary star 
is chosen randomly from a lKroupal (|2002l ) IMF of the form 



N(M) <x 



M 
M 



-1.3 
-2.3 



mo < M/M < mi . 
mi < M/Mq < m 2 . 



(X) 



where m = 0.1 M Q , mi = 0.5 M , and m 2 = 50 M Q . We 
do not include brown dwarfs in the simulations; the binary 
properties of brown dwarfs and very low mass stars appear 
to b e very different from t hose of M-, K-, and G -dwarfs 
(e.g. iBurgasser et~ai1 120071 ; iThies fc Kroupal 120081 ). For a 
fuller discussion of the effects of dynamical processing on 
brown dwarfs we re fer th e interested reader t o the work of 
iKroupa et ail (|2003l ) and lParker fc Goodwir] (|201ll ). 

We then assign a secondary component to the sys- 
tem depending on the binary fraction associated with the 
primary mass. For a field-like binary fraction we divide 
primaries into four groups. Primary masses in the range 
0.1 < M/M p < 0.47 are M-dwa rfs, with a binary frac- 
tion of 0.42 ((Fischer fc Marcvlll992l) . K-dwarfs have masses 
in the range 0.47 ^ M/M a < 0.84 with a binary frac- 
tion of 0.45 (|Mavor et al.lll992l ). and G-dwarfs have masses 
from 0. 84 ^ M/M B < 1.2 with a binary fractio n 
of 0.57 l|Duauennov fc Mavorl [l99ll ; iRaghavan etal1l2010l ). 
All stars more massive than 1.2 Mq are grouped together 
and assigned a binary fraction of unity, as massive stars 
have a much larger binary fraction than low-mass stars 
fe.g.lAbt et al.lll990l;lMason et al.lll998l;lkouwenhoven et all 
l2005l.l2007i ; |Pfalzner fc Olczakll2007l ; iMason et al.ll2009l . and 
references therein). 

We also set up clusters with a binary fraction of unity 
for all stars, and a binary fraction of 0.75 for all stars, ac- 
cording to the hypothesis that most, if not all stars, form in 
binary systems and that single stars are purely the result of 
dynamical processing of binaries and hig her-order systems 
l|Kroupalll995al ; IGoodwin fc Kroupall2005 
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2.2 Binary properties 

Secondary masses a re drawn from a fla t mass ratio distribu- 
tion; recent work by Reg giani fc Meverl (|201l[ ) has shown the 
companion mass ratio of field stars to be consistent with be- 
ing drawn from a flat distribution, rather than random pair- 
ing from the IMF. Currently, however, there is no detailed 
statistical analysis for the ONC. We note that drawing com- 
panions from a flat distribution means we do not recover a 
Kroupa IMF. 

We draw the periods of the binary systems from two 
generating functions. Firstly, in accordance with observa- 
tions of the fie ld, we use the logio-normal fi t to the G-dwarfs 
in the fiel d bvlDuquennov fc Mayorl l|l99ll . hereafter DM91) 
- see also iRaghavan et al.l l|2010l ). which has also been ex- 
trapola ted to fit the period distributions of the K - and M- 
dwarfs jMavor et al.lll992! ; lFischer fc Marcvll 19921 ): 



/ (logioP) oc exp 



-(logioP - logio-P) 

2(T 2 



(2) 



where logioP = 4.8, o~i oglo p = 2.3 and P is in days. Alterna- 
tively, we draw periods fr om the initial p re-main sequence 
period function derived by Kroupa ( 1995a. b, hereafter K95): 



/ (log 10 P) = V 



lQglpP - loglpPmin 
5 + (logioP - logloPnin 



(3) 



where logioPnm is the logarithm of the minimum period 
in days. We adopt logioPmm = 0; and r? = 3.5 and S — 
100 are the numerical constants ado pted bv lKroupal (|l995al ) 
and iKroupa fc Petr-Gotzensl (|201ll ) to fit the observed pre- 
main sequence distributions. We convert the periods to semi- 
major axes using the masses of the binary components. 

The eccentri cities of binary stars are draw n from a ther- 
mal distribution l|Heggielll975l ; iKroupal 120081 ) of the form 

Me) = 2e. (4) 

In the sample of iDuauennov fc Mayor! l|l99lf) . close binaries 
(with periods less than 10 days) are almost exclusively on 
tidally circularised orbits. We account for this by reselecting 
the eccentricity of a system if it exceeds the following period- 
dependent valucQ: 



e tid = 77 [0.95 + tanh (0.6 logioP - 1.7)] 



(5) 



We combine the primary and secondary masses of the 
binaries with their semi-major axes and eccentricities to de- 
termine the relative velocity and radial components of the 
stars in each system. The binaries are then placed at the 
centre of mass and velocity for each system in the fractal. 
The simulations arc run for 10 Myr using the kira integra- 
tor in the Starlab package (e.g. iPortegies Zwart et al.|[l999l . 
2001). We do not include stellar evolution in the simulations. 
Details of each simulation are presented in Table [T] 



1 iKroupal l ll995bl) and lKroupal teOOSft provides a more elaborate 
'eigcnevolution' mechanism to incorporate interactions between 
the primary star and its protostellar disk during tidal circulari- 
sation. However, this mechanism also alters the mass ratio distri- 
bution, causing a deviation fr om the flat mass ratio d istribution 
observed in the Galactic field jReggiani fc MevedliuTlT) . 



Table 1. A summary of the different cluster properties adopted 
for the simulations. The values in the columns are: the number 
of stars in each cluster (Altars), the typical mass of this cluster 
(M c i ustcr ), the initial fractal dimension of the cluster (D), the 
initial binary fraction in the cluster (/bin); aEL d the primordial bi- 
nary sepa ration di s tributi on (either iDuqucnno v fc Mavorl lll99ll . 
DM91) or lKroupal Jl995al . K95)). 



■^V stars 


-^cluster 


D 


/bin 


Separations 


1500 


~ 10 3 M 


1.6 


100 per cent 


DM91 


1500 


~ 10 3 M 


1.6 


100 per cent 


K95 


1500 


- 10 3 Mq 


1.6 


field-like 


DM91 


1500 


~ 10 3 Mq 


1.6 


75 per cent 


DM91 


1500 


~ 10 3 M 


2.0 


100 per cent 


DM91 


1500 


~ 10 3 Mq 


2.0 


100 per cent 


K95 


1500 


~ 10 3 Mq 


2.0 


field-like 


DM91 


1500 


~ 10 3 Mq 


2.0 


75 per cent 


DM91 


1500 


~ 10 3 Mq 


3.0 


100 per cent 


DM91 


1500 


~ 10 3 Mq 


3.0 


100 per cent 


K95 


1500 


~ 10 3 Mq 


3.0 


field-like 


DM91 


1500 


~ 10 3 Mq 


3.0 


75 per cent 


DM91 



3 RESULTS 

In this section we will describe the results of dynamical evo- 
lution on the primordial binary population in subvirial clus- 
ters with three differing levels of substructure, as set by the 
fractal dimension. We consider clusters with an initial frac- 
tal dimensions of D = 1.6 (highly substructured) D — 2.0, 
and D = 3.0 (almost no initial substructure). We first exam- 
ine the evolution of the substructure in the clusters, before 
following the evolution of the binary populations by looking 
at the overall binary fractions and separation distributions. 

We determine whether a star is in a bound binary 
system using the near est-neighbour method out lined in 
iParker et all (|2009h and iKouwenhoven et all i|2010l ). 



3.1 Cluster morphologies and evolution 

In Fig. [T] we show typical examples of initial cluster mor- 
ph ologies for the t hree initial levels of substructure. As found 
by I Allison et al.l (|2009l . [20101 ). the clusters collapse on very 
short timescales (< lMyr), leading to Plummer sphere-like 
morph ologies on tim escales of order the age of the ONC 
( 1 Mvr: I Jeffrie J2007al ibh . Irrespective of the initial fractal di- 
mension, the clusters reach similar morphologies after 1 Myr. 

In Fig. [2] we show the evolution of a typical cluster 
'coreQ over the lifetime of the simulation. We have picked 
a D — 2.0 simulation, but clusters with different fractal 
dimensions exhibit very similar behaviour. In this figure we 
plot the core density of the cluster as a function of time. The 
initial density of the fractal is 330 Mq pc -3 , which increases 
to 1920 Mq pc -3 during the densest phase at 0.9 Myr, im- 
mediately after cool-collapse. Following this dense phase the 



2 Note that it is meaningless to define a 'core' for a fractal cluster 
before it undergoes collapse. Before the collapse and formation of 
a core, we calculate the density within the half-mass radius from 
the centre of mass of the cluster. 



4 R. J. Parker, S. P. Goodwin and R. J. Allison 



1 - 



o - 



-2 



-2 





x/pc 



- 



-2 



-2 



: '. ••' i -Vj*W',-V 





x/pc 



(a) D = 1.6, Myr 



(b) D = 1.6, 1 Myr 



'■■■Tit.*- 



'fi ■mi. ' 



>v •••• 





x/pc 



o 

x/pc 



(c) D = 2.0, Myr 



(d) D = 2.0, 1 Myr 





x/pc 



(e) D = 3.0, Myr 



■-■••mam,- .. •• 



.■•.•-■*SS£»'*-» 
..'vKii.*-:- .: 



o 

x/pc 



(f) D = 3.0, 1 Myr 



Figure 1. Morphologies of typical examples of our clusters with an initial fractal dimension D = 1.6 at (a) OMyr and (b) 1 Myr; D = 2.0 
at (c) OMyr and (d) 1 Myr; and D = 3.0 at (c) OMyr and (f) 1 Myr. 
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Figure 2. The core density of a cluster with fractal dimension 
D = 2.0, undergoing cool collapse, as a function of time. The 
cluster reaches a peak density of 1920 Mq pc — 3 at 0.9 Myr. 



cluster quickly relaxes and after 10 Myr has a density of 
25 M pc~ 3 . 



3.2 Evolution of the binary fraction 

In Fig. |3] we show the evolution of the binary fraction over 
10 Myr, averaging together 10 clusters with the same ini- 
tial fractal dimensions. We show the evolution of the binary 
fraction for four different primordial binary populations; the 
DM91 distribution with an initially 100 per cent binary frac- 
tion (the solid line) , the K95 separation distribution with an 
initially 100 per cent binary fraction (the dashed line), and 
the DM91 separation distributions with field-like and 75 per 
cent initial binary fractions (the dot-dashed and dotted lines, 
respectively). 

The results shown in Fig. [3] are summarised in Ta- 
ble [2] For the various initial conditions, we show the bi- 
nary fraction as measured by our algorithm at OMyr, 1 Myr 
and 10 Myr. When considering the evolution of the binary 
fractio n in dense, virialised Plummer spheres, IParker et al.l 
(2009) noted that the cluster was too dense initially for the 
widest binaries observed in the field to be bound systems. 
T his means tha t the initial binary fraction in the clusters 
m IParker et af] {2009) with a 100 per cent primordial bi- 
nary fraction actually translated into an initial value of 75 
per cent for a DM91 separation distribution. The fractal 
clusters presented here are less dense than these Plummer 
spheres initially, and the calculated binary fractions are all 
higher (although none are 100 per cent). 

For clusters with a moderate level of substructure (D — 
2.0), an initial input binary fraction of 100 per cent, and sep- 
arations drawn from the DM91 distribution, the measured 
binary fraction at OMyr is 83 per cent, higher than the 75 
per cent initial bin ary fraction in a dense Plummer sphere 
l|Parker et al.ll2009T ). 

This effect is even more pronounced for the clusters 
with a binary fraction of 100 per cent and separations drawn 
from the K95 distribution. This distribution was derived to 
reconcile the observed overabundance of wide binaries in 
young clusters with the DM91 field distribution. Recently, 



Table 2. A summary of the results presented in Fig. [3] From left 
to right, the fractal dimension of the cluster (D), binary separa- 
tion distribution (/ (logioP)), initial binary fraction inputted into 
the simulations (/bin i n it); the initial binary fraction as measured 
by our algorithm (/bin,0Myr)i the binary fraction after 1 Myr 
(/bin,l Myr)i ar >d the binary fraction after 10 Myr (/bin, 10 Myr)- 



D 


/(logioP) 


/bin, init 


/bin, Myr 


/bin, 1 Myr 


/bin, 10 Myr 


1.6 


DM91 


1.00 


0.79 


0.55 


0.51 


1.6 


K95 


1.00 


0.68 


0.40 


0.38 


1.6 


DM91 


0.75 


0.62 


0.44 


0.40 


1.6 


DM91 


0.45 


0.41 


0.30 


0.28 


2.0 


DM91 


1.00 


0.83 


0.58 


0.54 


2.0 


K95 


1.00 


0.78 


0.48 


0.43 


2.0 


DM91 


0.75 


0.66 


0.51 


0.46 


2.0 


DM91 


0.45 


0.42 


0.34 


0.31 


3.0 


DM91 


1.00 


0.89 


0.74 


0.66 


3.0 


K95 


1.00 


0.88 


0.58 


0.51 


3.0 


DM91 


0.75 


0.68 


0.57 


0.53 


3.0 


DM91 


0.45 


0.44 


0.38 


0.35 



iMarks. Kroupa fe Obi (|201 lh have suggested that a dynam- 
ical operator (which is a function of the cluster's density) 
can be used to transform a K95 distribution to the field dis- 
tribution in a dense cluster. However, the K95 distribution 
saturates a dense cluster with wide binaries which are not 
physically bound, and it is dif ficult to see h ow th ey could 
form in such an environment. IMarks et alJ l)201ll ) suggest 
this problem could be negated if the cluster formed in a 
more sparse environment, and then underwent cool collapse, 
which is exactly the scenario we propose here. However, the 
calculated initial binary fraction in all the clusters here is 
significantly lower than 100 per cent (the dashed lines in 
Fig. [3] see also Table [2]| , which indicates that very wide bi- 
naries cannot form in star forming regions; an alternative 
solution is that they form during cluster dissolution, when 
two stars are simultaneously ejected in the same dir ection 
(e.g. iKouwenhoven et al. I l20ld : " iMoeckel fc Batdlioioh . 

The initial and final binary fractions depend heavily on 
the level of substructure. Comparing the simulations with 
D — 1.6 (highly substructured), to those with D = 3.0 
(uniform spheres), we see that the initial binary fraction 
is higher by 10 per cent for the uniform sphere (0.89 versus 
0.79), and after 10 Myr the difference is still significant, with 
a binary fraction of 0.66 for the D = 3.0 model versus 0.51 
for D = 1.6 model. 



Indeed, comparison of Figs. 3(a) and 3(b) with the over- 
all evolution of the cluster in Fig. [2] shows that the vast ma- 
jority of binary processing occurs before the cluster reaches 
its densest phase (after 0.9 Myr). This is due to pockets of lo- 
calised density in the substructure, which dynamically pro- 
cess the binary populations. In the case of a cluster with 
almost no initial substructure (Fig. 3(c) I, we see that there 
is very little binary processing until the cluster has almost 
reached its densest phase at collapse (note the sudden drop 
in binary fraction between 0.3 and 0.9 Myr, which corre- 
sponds to the density peak in Fig. [2} . 

The fact that dense substructure processes binaries to 
almost the same extent as the overall collapse of the clus- 
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Figure 3. The evolution of the binary fraction in clusters with different amounts of substructure; (a) a very clumpy cluster (fractal 
dimension D = 1.6), (b) a moderately substructured cluster (D = 2.0) and (c) a roughly uniform sphere (D = 3.0). Four different 
primordial binary populations are shown; (i) an initially 100 per cent binary fraction with the DM91 separation distribution (the solid 
line), (ii) an initially 100 per cent binary fraction with the K95 pre-main sequence separation distribution (the dashed line), (iii) a DM91 
separation distribution with an initially field-like binary fraction (the dot-dashed line), and (iv) a DM91 separation distribution with an 
initially 75 per cent binary fraction (the dotted line). 



ter suggests that substructured clusters in virial equilibrium 
and those undergoing expansion would also process any pri- 
mor dial binary popul ation. 

iPetr et all (|l998T l and iReipurth etHl l|2007t ) estimate 
that the binary fraction in the ONC is consistent with the 
field value, i.e. between 40 and 60 per cent, depending on 
the spectral type of the primary. We note from Fig. [3] and 
Table[2]that dynamical processing reduces the overall binary 
fraction to such an extent that the initial binary fraction 
cannot be that of the field. Even in the smooth clusters, the 
binary fraction at 1 Myr is less than 40 per cent. For the 
other initial conditions (D = 1.6 or D = 2.0), the binary 
fraction in the ONC at 1 Myr can be reproduced (within 
the uncertainties) if the initial binary fraction was 75 per 
cent or higher. 



ditions, a significant amount of dynamical processing re- 
duces the number of intermediate binaries, leading to an 
ov erall deficit of system s compared to the field. As found 
by iKroupa et alj (1 19991 ). the K95 separation distribution is 
reduced by interactions to the extent that the resultant sep- 
aration distribution resembles that of the field for close and 

intermediate separation bi naries. 

However, a s noted by iParker et af] 1 20091). and subse- 
quen t authors jKouwenhoven et al. 201C ; Moeckel fc Bate! 
bold ), no cluster undergoing a dense phase will pre- 
serve the wide binary systems observed in the field, 
and ot her mechanisms are requi re d to explain such sys- 
tems jKouwenhoyen et alj l20ld : IMoeckel fc Bate! l20ld ; 
iMoeckel fc Clarkejl201ll ~ 



3.3 The complete binary separation distribution 

We evolve clusters with two different initial separation dis- 
tributions. We consider clusters with separations drawn 
from the logio-normal distrib ution observed for main se- 
quence binaries in t he field (|Duquennov fc Mayor! Il99ll ; 
iRaehavan et alJl20ld ). and als o the in ferred pre-main se- 
quence distribution in iKroupal l|l995al ). Three out of four 
clusters have the DM91 separation distribution with vary- 
ing primordial binary fractions; 100 per cent, 75 per cent and 
field-like; whereas the final cluster has the K95 separation 
distribution. The initial separation distributions (the open 
histograms), and the distributions after 1 Myr (the hashed 
histograms) are shown in Fig. [4] For comparison we show the 
logio-normal fits to the separation distributions of field G- 
dwarfs (the (red) solid line), and field M-dwarfs (the (blue) 
dashed line). 

From inspection, we see that the results are simi- 
lar to those obtained with the virialised. dense Plummer 
sphere models presented in IParker et all (2009, see their 
figs. 2 and 3). In the model in which we use the field sepa- 
ration distribution and binary fraction as our initial con- 



3.4 Visual binaries in the ONC 

Furthermore, we note that the field binary population is 
proba bly the sum of many d i ffering star for mation re- 
gions (IBrandner fc Kohlerlll998l ; lGoodwin|[2oToh . not all of 
which will have undergone the cool-collapse scenario pre- 
sented here. Additionally, the most recent and complete 
observational census of binary systems in the ONC by 
IReipurth et al.l l|2007h only considers visual binaries with 
separations in the range 67 - 670 au. For this reason, it 
makes more sense to compare the results of our simulations 
to the data from Reipurth et al., rat her than the fie l d sep - 
aration distribution. The data from IReipurth et al l (|200Th 
are shown by the (green) crosses in Figs. [5] [6] and [71 We 
show our sepa r ation distributions (in the same range as 
IReipurth etall |2007l )1 after 1 Myr of dynamical evolution 
with the histograms (and the corresponding error bars from 
averaging together 10 simulations for each plot). 

In Figs.[5][6]and[7] we show the effects of dynamical evo- 
lution on the visual binaries in our clusters for the four initial 
binary populations. Firstly, we note that an initially field- 
like population (panel (c) in the figures) underproduces the 
required number of binaries in this separation range, apart 
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Figure 4. The full separation distributions at Myr (open histograms) and 1 Myr (hashed histograms) for clusters with an initial 
fractal dimension D = 2.0. Four different primordial binary population set-ups are presented, and the logio-normal fits to the separation 
distributions for G-dwarfs and M-dwarfs in the field (the (red) solid and (blue) dashed lines, respectively) are shown for comparison. 



from the clusters with smooth initial conditions (Fig. 7(c) I 



However, this cannot be the primordial binary population of 
the ONC because the overall binary fraction is lower than 



is observed (the dot-dashed line in Fig. 3(c) I. 

Secondly, and following on from this, the other separa- 
tion distributions for binaries in initially smooth clusters 
show that the number of visual binaries is overproduced 
(Figs. 



7(a) 7(b) and 7(d) I 



Finally, we see from inspection of Figs. [5] and [6] that 
all populations with an initial binary fraction of either 75 
per cent or 100 per cent reproduce the observed separa- 
tion distribution within the uncertainties, suggesting that 
there must have been an overabundance of binaries with 
separations in this range at the birth of the cluster. Be- 
cause of the highly uncertain binary fraction in the ONC 



|Petr et al J 1 19981 ; iKaczmarek et al1l201l[ ), we see from in- 
spection of Fig. [3] that clusters with an initial fractal di- 
mension of D — 2.0 or 1.6 are equally consistent with the 
observations, assuming either a DM91 or K95 initial sepa- 
ration distribution, and a binary fraction between 0.75 and 
unity. 



4 DISCUSSION 

We have examined the dynamical evolution of fractal clus- 
ters in cool collapse with three different levels of initial sub- 
structure. We consider clusters with fractal dimensions of 
1.6 (very clumpy), 2.0 and 3.0 (a roughly uniform sphere). 
In each of these substructured clusters, we examine the ef- 
fects of this cool collapse on four different primordial bi- 
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Figure 5. Comparison with the data for visual binaries in the ONC from Rcipurth ct al. (2007). Reipurth et al.'s data are shown 
by the (green) crosses with corresponding error bars. The separation distribution, normalised to the binary fraction at 1 Myr, from 
our simulations in clusters with an initial fractal dimension D = 1.6, are shown by the histograms. Four different primordial binary 
population set-ups are presented, and the logio-normal fits to the separation distributions for G-dwarfs and M-dwarfs in the field (the 
(red) solid and (blue) dashed lines, respectively) are shown for comparison. 



nary populations, which are characterised by the primordial 
binary fraction and binary separation distribution. Obser- 
vations indicate that many young star forming regions are 
both substr uctured, and that stars form with subvir ial ve- 
locities (e.g. iPeretto et al1l2006l ; iProszkow et al. || 2009l). Star 
clusters with these characteristics were used bv lAllison et al.l 
to show that mass segregation in the ONC can 
occur dynamically on a very short timescale ( 1 Myr) , negat- 
ing the need for primor dial mass segr egation in the ONC 
(|Bonnell fc Davieg [T998). Furthermore, lAllison fc Goodwin! 



|201ll ) have shown that trapezium-like systems regularly 
form in such simulations, suggesting the cool-collapse of a 
clumpy cluster could be the most likely dynamical evolution 
scenario for the ONC. The most favourable initial conditions 
for this dynamical mass segregation (and the formation of 
the Trapeziu m system) are a clumpy (D 2.0), cool (Q < 
0.4) cluster (|Allison et al ] 120091 , 120101 : lAllison fc Goodwin! 

The hypothesis presented in lAllison et al.l l|2009l ) is sup- 
ported by observations. The outskirts of the ONC (20 pc 
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(c) Field binary fraction, DM91 separation distribution 



(d) 75 per cent binary fraction, DM91 separation distribution 



Figure 6. Comparison with the data for visual binaries in the ONC from Reipurth ct al. (2007). Reipurth et al.'s data are shown 
by the (green) crosses with corresponding error bars. The separation distribution, normalised to the binary fraction at 1 Myr, from 
our simulations in clusters with an initial fractal dimension D = 2.0, are shown by the histograms. Four different primordial binary 
population set-ups are presented, and the logio-normal fits to the separation distributions for G-dwarfs and M-dwarfs in the field (the 
(red) solid and (blue) dashed lines, respectively) are shown for comparison. 



from t he Trapezium) appear to be subvirial a nd in cool col- 
lapse |Feigelson et al.ll2005l ; iTobin et al.ll2009h . whereas the 
velocity dispersion in the centre is 4.3 kms -1 , much higher 
than the value we wo uld expect if the ON C was in virial equi- 
librium (2.5 kms" 1 - lOlczak et~al] [2008). This suggests that 
the centre of the ONC has already undergone cool-collapse, 
and is now expanding. 



However, for simplicity I Allison et al.l (120091 ) did not in- 
clude primordial binaries in their simulations. The binary 
fraction in the ONC is not negligible, and is consistent with 



the field value (betw een 40 and 60 per c entlPetr et al] 



Reipurth et al. 2007 ). Several authors jKroupa et al 



1998 



1999; 



Parker et alj 20091 ; iKaczmarek et alj l201ll ) have proposed 
that the ONC was born with a much higher binary frac- 
tion than its present value, and that dynamical interactions 
have processed this primordial population to that which we 
observe today. For this to happen, the cluster must have un- 
dergone a dense phase during its evolution. If the results of 
dynamical processing on a primordial binary population can 
be reconciled with the observed binary fraction and separa- 
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Figure 7. Comparison with the data for visual binaries in the ONC from Rcipurth ct al. (2007). Reipurth et al.'s data are shown 
by the (green) crosses with corresponding error bars. The separation distribution, normalised to the binary fraction at 1 Myr, from 
our simulations in clusters with an initial fractal dimension D = 3.0, are shown by the histograms. Four different primordial binary 
population set-ups are presented, and the logio-normal fits to the separation distributions for G-dwarfs and M-dwarfs in the field (the 
(red) solid and (blue) dashed lines, respectively) are shown for comparison. 



tion distribution, then this provides strong support for such 
a theory. However, it is unclear how a significant nu mber 
of binaries could form in such a dense environ ment (|Batei 
120091 ; iParker et all [20091 ; iMoeckel fc Bate! l201Ch . An alter- 
native scenario is that the initial density of the cluster is 
such that binary formation is not impeded, but the cluster 
then undergo es a dense phase via the c ollapse of a substruc- 
tured fractal (| Allison et al Jl200Sl . l201oh ■ An excellent test of 
this hypothesis is to study the effects of cool-collapse on a 
primordial binary population. 



The low, field-like binary fraction in Orion is in itself not 
conclusive proof that the cluster has undergone significant 
dynamical processing; a more stringe nt test is to exam i ne the 
separation distribution of the cluster. iReipurth et~ai1 (|2007f ) 
conducted a survey of visual binaries in the ONC, corre- 
sponding to a separation range 67 - 670 au. If the ONC did 
go through a dense phase, then we would expect the hard - 
soft boundary for binary disruption jHeggiell975l ; lHillslll975l ) 
to lie within this range l|Parker et al.ll2009l ). Therefore, by 
examining the effects of cluster evolution on various primor- 
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dial binary populations, we can constrain the primordial bi- 
nary fraction and separation distributi on (in this separation 
range) based on comparison with the lReipurth et aD (|2007l ) 
data. 

Direct comparison of the observations with our simu- 
lations is presented for each initial level of subclustering 
in Figs. [5l [6] and [7] Firstly, we note that the process of 
cool-collapse in clusters can reproduce the observed separa- 
tion distribution to zeroth order for most primordial binary 
populations. Clusters with a moderate to high level of sub- 
clustering cannot preserve enough binaries in the separation 
range 67 - 670 au for an initially field-like binary fraction 



(Figs. 6(c) and 5(c) I. If we start the cluster as a uniform 



sphere (D = 3.0), then it is possible to reproduce the obser- 
vations with a field-like binary fraction and separation dis- 



tribution (Fig. 7(c) I. However, the overall binary fraction for 
the cluster is still too low (Fig. 3(c) I, suggesting that even in 
this more placid dynamical scenario, the primordial binary 
fraction has to be larger than the present day. Furthermore, 
all other initial binary populations do not undergo enough 
processing to suggest that this fractal dimension is a realistic 
initial condition for the ONC. 

A moderate level of of substructure (D — 2.0) results in 
excell ent agreement with the observations of lReipurth et all 
for clusters with DM9f separation distributions 
and primordial binary fractions of f00 or 75 per cent 
(Figs. |6(a)| and |6(d)| respectively). In clusters with very 
clumpy initial conditions (D — 1.6), the level of dynamical 
processing is too extreme in all but the cluster with a DM91 
separation distribution and a 100 per cent primordial binary 
fraction (Fig. 5(a)| to be reconciled with the observations of 
iReipurth et ;uTp007h . 

As discussed in lAllison et all ( |2010h and 
lAllison fc Goodwinl (|201 lh . dynamical mass segrega- 
tion and the formation of trapezium systems can be very 
transient. In order to reproduce the observed level of 
mass segregation it is favourable to have clumpy initial 
conditions. If we assume that clumpy, cool initial conditions 
are required for the ONC to mass segregate and form the 
Trapezium system, then the observed binary fraction and 
separation distribution requires an initially higher binary 
fraction (~ 70 - 80 per cent) than is observed today. 

Finally, we note that if the clusters are initially clumpy, 
the majority of binaries are processed before the cluster 
reaches its densest phase during the collapse. This is be- 
cause the pockets of substructure are dense enough initially 
to affect the binaries, and suggests that all star clusters that 
form with substructure will process a primordial binary pop- 
ulation, irrespective of whether the cluster undergoes cool- 
collapse (which exacerbates the processing) , remains in virial 
equilibrium, or expands. 

Therefore, even some expanding associations which 
form supervirial/unbound may not preserve their primor- 
dial binary populations. However, we note that the frac- 
tals we set up have initial densities of ~ 300 Mqpc -3 (see 
Fig- HJ, which are higher than many star forming regions 
that will subsequent l y become unbound as s ociations (e.g. 
Jorae nsen et al . 2008; iGutermuth et al . 20091; iBressert et al.l 



20101 ). We will further investigate the effects of dynamical 
evolution on such sparse regions in a future paper 



5 CONCLUSIONS 

We present the results of A^-body simulations of fractal star 
clusters containing N = 1500 stars in cool-collapse, in order 
to investigate the effect of this dynamical evolution scenario 
on various primordial binary populations. We have varied 
the initial level of substructure in the cluster, the primordial 
binary fraction, and the initial separation distribution. Our 
conclusions can be summarised as follows: 

(i) Primordial binary populations are heavily processed 
in clusters undergoing cool-collapse. Qualitatively, the re- 
sults are similar to those from dynamical evolution of the 
binary p opulation in initially very dense virialis ed Plummer 
spheres jKroupa et al.lll999l ; |Parker et alj|2009h . 

(ii) The level of dynamical processing varies as a func- 
tion of the fractal dimension; clumpy clusters break up more 
binaries than smoother clusters. 

(iii) The majority of dynamical processing in substruc- 
tured clusters occurs before the cluster reaches its densest 
phase; therefore, it is the initial densities in the substruc- 
ture which is the most significant contributor to altering the 
binary population, rather than the cool-collapse itself. This 
suggests that even some star-forming regions that do not 
collapse will significantly process a primordial binary popu- 
lation. 

(iv) If clusters undergo cool-collapse, then the field bi- 
nary fraction and separation distribution cannot be the pri- 
mordial distribution in the ONC. Comparison of our simu- 
lations with observations suggests that the ONC had a pri- 
mordial binary fraction of between 75 and 100 per cent. 

We demonstrate that the cool-collapse scenario, which 
is consistent with the filamentary, subvirial early phases of 
star formation, and can explain the level of mass segregation 
in the ONC through dynamics, also reproduces the observed 
binary fraction and separation distribution. If a moderate to 
high level of substructure is required to produce dynamical 
mass segregation, then an ~80 per cent binary fraction, and 
field-like separation distribution with a cut-off around 5 x 
10 3 au represents the most likely initial binary population. 
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